CMB anisotropies in the presence of a stochastic magnetic field 
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Abstract 



Primordial magnetic fields present since before the epoch of matter-radiation equality have an 
effect on the anisotropies of the cosmic microwave background. The CMB anisotropies due to 
scalar perturbations are calculated in the gauge invariant formalism for magnetized adiabatic initial 
conditions. Furthermore the linear matter power spectrum is calculated. Numerical solutions are 
complemented by a qualitative analysis. 



1 Introduction 

> 

£C) ' Observations of the cosmic microwave background (CMB) provide an important tool for testing 

^& the understanding of the physics of the early universe. Since the first data of the angular power 

spectrum of the temperature fluctuations in the CMB provided by COBE [I] the aim has been 
to put limits on cosmological parameters such as its matter composition or present day expansion 
rate. Moreover, with more and better quality data not only of the temperature fluctuations but 
also the polarization of the CMB becoming available at a fast rate it is possible to test even earlier 
stages of the universe long before the beginning of the standard radiation dominated era. Thus 
different models of inflation are under scrutiny and part of the class of models has already been 
S^ . found to be excluded by observations [2] . 

Apart from matter, radiation and dark energy there is evidence for magnetic fields on small 
upto very large scales. Magnetic fields on galactic scales have a field strength of the order of a /iG 
and correlation lengths of the order of few kpc. There are also observations of magnetic fields in 
clusters which have similar field strengths and correlation lengths. More recently there have been 
indications of the existence of magnetic fields in high redshift galaxies [3 J • 

There exist a multitude of proposals of generation mechanisms of magnetic fields on scales 
comparable to galactic scales. Generally speaking there are two types of generation mechanisms. 
Firstly models where magnetic fields are generated in the very early universe during inflation due 
to the amplification of perturbations in the electromagnetic field. Secondly battery-type models 
operating after inflation has finished [3]. 

In most models magnetic fields are present long before the initial conditions for the evolution 
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of the metric and matter perturbations are set, significant for the formation of the characteristic 
spectrum of the anisotropies of the CMB and its polarization. Therefore taking into account 
the presence of such a primordial magnetic field provides an interesting tool to put limits on 
its characteristics such as its field strength or spectral index. This has been already studied in 
different settings. The effect of the primordial magnetic field on the anisotropies of the cosmic 
microwave background (CMB) depends on its nature. Assuming the magnetic field to be stochastic 
guaranties that it does not break the global isotropy of space-time [4]. Moreover there have also 
been investigations where the primordial magnetic field is assumed to be uniform, which for large 
values of the field strength would lead to models incompatible with the observed large scale isotropy 
of the universe [5] . Here we are going to assume that the magnetic field is primordial in origin and 
thus present since long before matter-radiation equality and, moreover, it is a stochastic magnetic 
field, so that global isotropy is preserved. 

There are several approaches to formulate the perturbation equations depending on whether a 
particular gauge is chosen, such as for example the synchronous gauge or a covariant formulation 
which is based on the fluid description. In this latter approach no particular gauge is chosen. The 
calculation of the anisotropies of the temperature and the polarization of the CMB requires to solve 
the Boltzmann hierarchy describing the evolution of the photon distribution function. COSMICS 
[6] was the first numerical Boltzmann solver which was followed by CMBFAST [7] which uses 
the synchronous gauge to describe the perturbation equations and introduced the line-of-sight 
integration speeding up significantly the numerical calculation of the CMB anisotropies. CAMB |8j 
is another numerical code to solve the Boltzmann hierarchy and determine the CMB anisotropies. 
It uses the covariant approach to formulate the perturbation equations. More recently CMBEASY 
[9] yet another Boltzmann solver has been released. It allows to calculate the CMB anisotropies 
using the gauge invariant formalism to formulate the perturbation equations. Apart from the 
different mathematical approaches used to formulate the perturbation equations it is interesting 
to note the evolution from the point of view of programming languages. Whereas COSMICS and 
CMBFAST are written in Fortran 77, CAMB uses Fortran 90/95 and CMBEASY is written as 
a completely object-oriented C+- 1- programme. The effect of a primordial magnetic field on the 
CMB anisotropies has been calculated using different approaches: synchronous gauge and thus a 
modified version of CMBFAST in |10t [TT] or the covariant formalism and a modified version of 
CAMB in pi [1311I4]. 

Here the gauge-invariant formulation of the perturbation equations is used and a modified 
version of CMBEASY to calculate the CMB anisotropies due to the scalar perturbations of the 
geometry. 

In section 2 the perturbation equations and the initial conditions in the gauge-invariant formal- 
ism are presented. Section 3 is devoted to the calculation of the different magnetic field contributions 
which due to the stochastic nature of the magnetic field involves convolution integrals. The mag- 
netic field spectrum is damped due to diffusion on small scales. This is implemented here by using 
a gaussian window function. In section 4 the angular power spectra determining the temperature 
and polarization auto-and cross-correlation functions and the linear power matter spectrum are 
presented. Section 5 contains the conclusions. 



2 Gauge-invariant description 

In the gauge-invariant formalism the perturbation equations are written in terms of gauge-invariant 
variables. For the scalar sector Einstein's equations in Fourier space are given by [15(. [TE] 

^pA + ^pjl + w^k-'V 

2M 2 k 2 + 3a 2 (l + w)p { ' ' 

* = -*- |§ '-) 

where $ and \1/ are the gauge-invariant Bardeen potentials [17] . A dot indicates the derivative 
with respect to conformal time of the unperturbed metric of the background space-time, that is 
ds 2 = a 2 (T)(—dT 2 + Sijdx l dx : >), and a(r) is the corresponding scale factor. A is the gauge invariant 
density perturbation defined by 

A = 6 + 3(l + w){Q+Hk- 1 a), (2.4) 

where 5 = ^^- is the coefficient of the density contrast in the harmonic expansion and a is the 
coefficient of the shear in the harmonic expansion. Photons (7), neutrinos (u), baryons (b), cold 
dark matter (c) and the magnetic field (B) contribute to the total energy density perturbation, so 
that 

pA = /) 7 (A 7 + A B ) + p i ,A I/ + p c A c + p b A b . (2.5) 

The total contribution to the gauge- invariant velocity V [151 [T6] and the total anisotropic stress 
are given by 

4 

(l+ W )pV = -(PyVy + PvVj+PcVc+pbVb, (2-6) 

P U = gP 7 ( 7r 7 + 7r B) + -PvK u . (2.7) 



Finally, 



4 
(1 + w)p= -{p~/ + Pu) + Pb + Pc- {2.1 



Furthermore Mp = Mp/y8ir is the reduced Planck mass and % = - leading to the form of the 
Friedmann equation, 



* 2 =sr (2 - 9) 

Well within the radiation dominated epoch electrons and baryons, coupled by Coulomb interaction, 
are as well tightly coupled to the photons because of Thomson scattering of the photons off the 
free electrons. Thus the baryon-electron-photon system is very well approximated by a one- fluid 
description. Due to the presence of the magnetic field this baryon-electron-photon fluid is magne- 
tized. As shown in |18j the system is in the magnetohydrodynamic limit. The baryons are subject 



to the Lorentz force which changes the evolution of the baryon velocity in Fourier space to (cf., 

e.g., usi nana nam) 



V h = (3c 2 - l)HV h + fc(* - 3c 2 s $) + kc 2 s A h + Rt" 1 ^ - V h 



+ -kL, 



(2.10) 



where L is due to the Lorentz force J x B and R = , — . Furthermore, ci 

sound speed and r" 1 is the mean free path of photons between scatterings given in terms of the 



25? is the adiabatic 



number density of free electrons n e and the Thomson cross section or, t c 



an e aT- In section 



[3] the contribution of the magnetic field given by its energy density, anisotropic stress and Lorentz 
term will be specified. 

2.1 The tight-coupling limit 

In the very early stages, long before recombination, the energy density of free electrons scales as 
a -3 so that the time between scatterings of photons is proportional to a 2 . Thus the mean free 
time between scatterings of the photons is much smaller than the Hubble time which leads to a 
comparatively large value of r" 1 . This implies that in the numerical integration of the baryon and 
the photon velocities special care has to be taken with respect to the time step size. This problem 
was solved by using an iterative solution during the very early stages in the tight coupling limit 
and afterwards the original equations |20(, [2T| I22j . To derive the equations for the evolution of 
the photon and the baryon velocities in the tight coupling limit we start by giving the relevant 
equations determining the photon evolution [151 IS] > 



v; 
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--kv, 

*(* - *) + ^A 7 - ^vr 7 +r~\V h - VJ 



(2.11) 
(2.12) 



and the baryon density contrast evolves as, 



A, 



-kV h - 3c 2 nA h . 



(2.13) 



The baryon velocity is determined by equation (|2.10p . From these equations using that R ~ a 
and c 2 ~ a -1 the so-called slip equation is derived V = Vh — V 7 [22] which is found to be, 



V 



1 + 2°- Tc 



al + R 



Tr 



A, 
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1 + R 

TT, 
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V h 



Tr. \ 2 6 



U ! c 2 A h - — + * + ^ 



+ -(v h - v 1 



(2.14) 



In the tight coupling limit the term Vb — V~ is neglected [22 . In terms of V the equations for the 
photon and baryon velocities are given by, 



v; 



R 



1 + R 
1 



-k 
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1 + R 



A~ 



f + 4 



$ + k^ 



H {3c 2 -l)V h + kc 2 (A b - 3$) - V 



(2.15) 



and 



Vb = — ^ l n ( 3c * ~ l ) v * + kc l ( A b - 3*)] + ^ 



1 + 22 

R 

+ 



l + R 



,,^-^ + |-<i>)+y 



(2.16) 



2.2 Initial conditions 



In the numerical calculation of the anisotropies in the CMB it is usual to set the initial conditions 
after neutrino decoupling so that the neutrino anisotropic stress is non zero. This led to the 
formulation of initial conditions for the calculation of the CMB anisotropies in the presence of 
a primordial magnetic field in which the contributions involving the neutrino anisotropic stress 
and the magnetic stress cancel each other at lowest order in x = kr |X8|, [TTj [T2j [T3"l [T4"] . These 
are the compensating initial conditions. However, in a large class of models the magnetic field 
is generated in the very early universe, long before neutrino decoupling, such as during inflation 
(e.g. |23j) or the electroweak phase transition (e.g. [24]). In this case the magnetic field provides 
initially the only source of anisotropic stress. It was shown in [25] for a general type of anisotropic 
stress and in |14|. [26] for the special case of a magnetic field that the solution for the neutrino 
anisotropic stress after neutrino decoupling approaches a solution compensating the contribution 
from the anisotropic stress of the magnetic field. Moreover, there is an additional contribution to 
the curvature perturbation on large scales. Therefore in [14] it was concluded that after neutrino 
decoupling there are effectively two types of perturbations associated with the primordial magnetic 
field. On the one hand there is the compensating mode and on the other hand there is the passive 
mode which is an adiabatic-like mode with a non-vanishing curvature amplitude. Therefore the 
approach that is taken here to find the numerical solutions is to set the initial conditions after 
neutrino decoupling and assume the compensating initial conditions. 

Following [16] new variables are defined by, Vi = Vi/x and 7fj = 7Tj/x 2 . Then the initial 
conditions on superhorizon scales, x <C 1, are determined by the set of first order differential 
equations, 

"' A --x 2 V **!L = -tfv ^--x 2 V ^L-- x 2 V h (2 17) 



dlnx 3 dlnx 3 ' dlnx dlnx 

dV^ A,v ~ „ „ _ L dV c 



"y 



dlnx 4 ^ + <¥fe+0,*, + 2* + I , — = -2F c + *(2.18) 

dV A x 2 

-p - V v + 2* + n 7 7t B + n„7T„ - -7T 7T„ (2.19) 



dlnx 4 '6 

dftv 8 ~ dnn 

-V u -2n u , ^L = _2# B . (2.20) 



dlnx 5 cilnx 

Solving this system to lowest order in x imposing A 7 = A u = |A C = |Ab, the following set of 



initial conditions is found, 



- 5 A 7 5 ft 7 (A B + L) 5 fi 7 3 - 2VL V 

v ~ ~415 + 4fi I/ ~2 15 + 40,, + 6f^l5 + 4^/ B 

5 A 7 5 7 A B 5 + 140,, L 7 7 tt b 



4 15 + 40^ 2 15 + 4fi„ 15 + 40^ 4 3 15 + A£l u 
V 5 A 7 5 - 4fi„ 7 13-40^0^ 

c 415 + 4fi„ 15 + 40,, 8 [ B J 15 + 4fi„ 12 

Q-y _ A 7 2Q 7 (A B + -L) 2 7 3 - 2Q. U 

*" ~T^ nB ~ 15 + 4fi„ 15 + 4^ + 3J^15 + 4n i/ 7rB " ^' ^ 

As shown in the appendix, these initial conditions correspond to the compensating initial conditions 
in the synchronous gauge used in previous numerical solutions |18|, lllj , fT2], IX3 |, I14j . Moreover, in the 
case of no magnetic field these (" magnetized adiabatic" ) initial conditions reduce to the standard 
adiabatic initial conditions |16j . 

The scalar curvature on a comoving hypersurface is given by (cf., e.g., |27j ) 

-£ = -^ + %k- l V. (2.22) 

The initial conditions are set after neutrino decoupling and we focus on the compensating magnetic 
mode. Therefore, using equation (I2.22p . initially the total comoving curvature perturbation is given 
by 

C = ^ + ^ (2.23) 

which can be used to express A 7 as A 7 = 4£ — Q 7 Ab- Moreover Q is treated as a Gaussian random 
variable with the two-point function in Fourier space (£*(&)£(£;')) = "P^S,^, an d the dimensionless 

n 2 f u \ n s — 1 
power spectrum is defined by Vq = "p-^4s I y~ I where n s is the scalar spectral index, A s the 

amplitude of the scalar perturbations and k p = 0.002 Mpc -1 the pivot wave number used in WMAP 
[28]. In the numerical solutions A s and n s will be set to the bestfit values of the six parameter 
ACDM model of WMAP7 [28]. Taking into account that the initial total curvature perturbation 
is given by the sum of the contribution resulting, for example, from inflation and the magnetic 
contribution, assuming £ to be determined by the bestfit values of WMAP and given the magnetic 
field parameters constrains the curvature perturbation from inflation. 

In \\.A\ [26] it was found that if the magnetic field is generated before neutrino decoupling the 
scalar curvature evolves on superhorizon scales to a final value at the time of neutrino decoupling 
determined by the magnetic anisotropic stress, the time of generation of the magnetic field tb and 
the time of neutrino decoupling r v . This leads to the passive magnetic mode which corresponds to 
an adiabatic-like mode with an amplitude given by 



S,pass — „^7^"B 



^khik- 1 



(2.24) 



where i? 7 = q t^ ■ If the magnetic field is generated during a phase transition, which is the case 
considered by [26J then t-q corresponds to the time of the phase transition, if however, it is generated 
during inflation t-q could be chosen to be at reheating [T5]. The passive mode has to be added to 
equation ()2.23j) which consequently modifies the constraint on the nonmagnetic contribution to the 
initial total curvature perturbation. 
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3 The magnetic field contribution 

In order to describe the magnetic field the so called "lab" frame is chosen, in which 

where eij m is the totally antisymmetric symbol, with £123 = 1 and F^ u is the Maxwell tensor. The 
lab frame is defined by choosing locally Minkowski space-time, so that the lab coordinates are 
defined by dt = adr, dr = adx [29]. Assuming that the conductivity of the plasma is large, Ohm's 
law implies a vanishing electric field, which leads to the magnetic field decaying as 1/a 2 as the 
universe expands. In general the energy momentum tensor of the electromagnetic field measured 
by the fundamental observer can be written in terms of that of an imperfect fluid [30] , 



T a p = (p + p)u a Up + pg afS + 2u (Qi g /3 ) + vr a/3 , (3.2) 

where u a = cl~ 1 5q is the 4-velocity of the fluid and u a u a = — 1. q a is the heat flux which, in 
the case of the electromagnetic field is determined by the Poynting vector and thus vanishes for 
vanishing electric field [30]. The magnetic energy density /?b, pressure pb and anisotropic stress 
7T(b) a/3 in the lab frame are given by [30] , 



B 2 (x,t) 1 1 

" 3PB, 7T( B )y =-B i (x,T)B j (x,T) + -, 



PB = 7T — > PB = oPB, ir {B)ij = -B i (x,T)B j (x,T) + -B (x,T)5ij, (3.3) 



where the anisotropic stress has only non-vanishing spatial components and the vector notation 
denotes a spatial 3-vector. Moreover, it has been used that the anisotropic stress tensor changes 
from the FRW frame to the lab frame as ttE 1 }^ = vri^-a 2 . Furthermore, the term due to the 
Lorentz force entering the equation of the baryon velocity evolution and in the tight-coupling limit 
the photon velocity evolution which is derived from v a T^ m = —F^ a J a and expressed in terms of 
quantities in the lab frame yields to 

L(x, r) = a(jxB) (x, r) (3.4) 



For vanishing electric field the current J is given by a J = V x B. Thus the components of the 
Lorentz term takes the form 

L J = -Q d iB 2 ~ Yl di7T (V) ij ( 3 - 5 ) 

i 

In the perturbation equations the magnetic field contributes to the total energy contrast and 
the total anisotropic stress. The magnetic energy density ps defines the magnetic energy density 
contrast 5b as 

Pb{x,t) = p^ J ^25 B (k)Y(k,x), (3.6) 

k 

where Y(k,x) denote a complete set of scalar harmonic functions satisfying, (A + k 2 )Y = |15j . 
The magnetic energy density is defined in such a way that it does not contribute to the total 
background energy density. This is different from [30] where the magnetic energy density does 
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contribute to the background. Moreover, Ab = #b- The anisotropic stress can be expanded in 
terms of scalar harmonics as, 

k 
where Yij = k~ 2 Y^j + ^6ijY [15]. Thus the Lorentz term can be written as 

Li(x, t) = ^Y. kL(k)Y(k, x) L(k) = A B - ^vr B (3.8) 

k 

and Yi = — k~ l Yu [15]. The magnetic energy density, the anisotropic stress and the Lorentz term in 
Fourier space can be related to the magnetic field spectrum. Using that B{(x, r) = Bi(x, To) ( ^y J 
and j0 7 = p y o (^-) , where the index refers to the present epoch and defining 

B i (x,T Q ) = Y,B i {k)Y{k,x), (3.9) 

k 

where Bi(k) = Bi(k,ro). Thus it is found that 

Ab(£) = ^~ V Biffi&ik - q) (3.10) 

<? 

and using the last equation of (|3.3p and equations f|3.T[) and (|3.9p yields 



7TB (k) 



2p 7 o 



Y, j^BMik' ~ j)Bj{% - q)qi - £ B m {q)B m {k 



(3.11) 



The magnetic field at present is characterized by its spectrum Pq which is chosen to be of the form 

PB{k,k m ,k L ) = A B ( r ) W{k,k m ), (3.12) 

where Ab is its amplitude, n^ is the magnetic spectral index and k m the upper cut-off in the 
magnetic field spectrum due to diffusion of the magnetic field energy density on small scales. &l is 
the pivot scale of the magnetic field and W(k, k m ) is the window function modeling the cut-off of 
the magnetic field spectrum due to diffusion. As shown in [29] the damping of the magnetic field 
is determined by the dimensionless Alfven velocity Va and the Silk damping scale kg, that is 

k m 2 = Vlk~ 2 , (3.13) 

where Va ~ 3.8 x W~ 4 (B/1 nG) and B corresponds to the smoothed field strength of the magnetic 
field today and before recombination the Silk scale may be approximated by [31] 



1.7 x 10' ( 1 - — ) (fi b k ) (fiok )~ 5 ( — — ; Mpc 2 , (3.14) 

. a oJ 3 v /a eq /a + 2 



where Yj> is the primordial helium mass fraction. After decoupling the photon mean free path 
becomes infinite and the viscous damping becomes subdominant. Therefore, the largest scale of 
the magnetic field spectrum damped corresponds to /c" 1 evaluated at recombination. Thus, it is 
interesting to note that the diffusion scale changes with time and this would lead to a magnetic 
field spectrum which has a more general time dependence than just the decay due to expansion. 
Therefore, since this case is not considered here, the magnetic damping scale is assumed to be 
defined by the largest damped scale. This scale is determined by the Alfven velocity and the Silk 
damping scale at recombination. The Silk damping scale is determined by the mean free path of 
the photons which continuously grows from the early epochs within the tight-coupling regime of 
the baryon-photon fluid to recombination shortly after which the photon mean free path becomes 
infinity. In |29j the maximal damped magnetic scale has been estimated and the corresponding 
maximal wave number is given by 

/c m ~ 200.694 ( — J Mpc" 1 , (3.15) 

for the values of the bestfit ACDM model of WMAP7, n b = 0.0227/i" 2 and h = 0.714 [28]. 
In the following the window function is assumed to be gaussian of the form, 

W(k,k m ) = ^-§£;- 3 eVW , (3.16) 

in such a way that J d 3 kW(k,k m ) = 1. This choice of window function is different from the step 
function used in previous work |1 1 (. [T2l [13"1 [T4"] . The smoothed magnetic field strength is defined 
by the auto-correlation function, (B 2 (x,tq)} which can be easily calculated using the correlation 
function of the magnetic field given by 

{B*{k)B 3 {q)) = S^Pb (<% - ^) , (3.17) 

where the helical part is neglected since for a stochastic magnetic field it does not contribute to 
the CMB anisotropies due to scalar perturbations [32]. Taking the continuum limit, Yl% ~~ ^ I (2-K)' i 
implies that the magnetic field strength today smoothed over the magnetic diffusion scale is given 
by 
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M nBr ( 2 ) 



{B(xY) = A**-* \jz.) -^y-A ( 3 - 18 ) 

which is valid for n-Q > —3. The damping scale is a natural scale in the problem. Therefore here 
it is not necessary to introduce a smoothing scale as it is usually done. Using equation f)3. 15j) the 
damping scale A m = 2tt /k m is of the order of A m ~ 30(-B/nG) kpc and thus in general smaller than 
1 Mpc which is usually used as smoothing scale. Defining the two-point correlation functions in 
/c-space in terms of the dimensionless spectrum Vf by 

2vr 2 

The following expressions are found for the spectra determining the auto correlation functions of the 
magnetic energy density contrast, the anisotropic stress and the Lorentz term. These are different 



from previously obtained expressions [T2J Q2J El E3] due to the different choice of window function. 
In general the spectrum characterizing the auto-correlation function of the magnetic energy density 
contrast is given by [I], 

2" 



k 3 



47r P-yO , 



1 



q-(k-q) 



Q 2 \k ~ Q\' 



(3.20) 



which for the spectral function (|3.12p leads to 

2 

VA B (k,k_ 



1 



[rmr 



PBO 



P7O 



2(n B +3) 



dzz nB 



+VH&J z 



ctee' 



[1 - 2zx + z 2 ) 2 (1 + :r 2 + 2z 2 - Azx) , (3.21) 



where x = k • q/(kq) and z = |. Moreover, the average energy density of the magnetic field is 
defined to be /?bo = {B 2 {x))/2. For small values of J*-, that is r-<l, the double integral can be 
approximated by incomplete Gamma functions. In particular, it is found that 



^A B (fc)^m) 



1 



2^ /* 



PBO 

P-yO 



-(riB+3) 



k 



2(n B +3) 



/ "B + 3 / fe 

\ 2 \ /Cm 



+ 2~( nB+ i) 



n'n 



-2n B — 3 



T U B + -;2 



n'li 



(3.22) 



where 7(0;, x) = J^ e~ t t a ~ 1 dt and T(a,x) = f°° e~ t t a ~ 1 dt are incomplete Gamma functions 
It is interesting to note that in the limit k <C k m the power-law behaviour which was used in 
the magnetic field spectra used in [TT] is recovered. However, as can be seen from figure Q] the 
approximation underestimates the contribution from the magnetic field and thus the effect on the 
final angular power spectrum of the temperature anisotropics and polarization. 

The spectrum defining the auto correlation function of the anisotropic stress is determined by 

12. 



?+ B (Mm) = -^^2P B (q)PB(\k-q\) 
87r P-yO a 



18q 2 {l-x 2 ) 2 q(kx - q)(l - x 2 ) 



k-q\ 2 

' (kx — q) 2 



\k-q\< 



+2 1 + 



q]' 



(3.23) 



Using the expression for the magnetic field spectrum (|3.12p and the window function (|3.16p this 
leads in the continuum limit to, 



' 7t b \ki kn 



dxe 



9 



(1 - 2zx + z 2 ) 



Pbo 
P70. 

"R-2 



k 



2(n B +3) 



dzz nB+2 e 2 W " 2 



1 + 5z 2 + 2zx + (1 - 12z 2 )x 2 - 6zx 3 + 9z 2 x 4 ) . (3.24) 
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This gives in the limit for J-<1, 



I n B \k, Km) 



[rm] 2 



PBO 

P70 



7±f" B+ Vte)" 

fern/ 



2~ 



n'll 



-(n B +3) 



H-b +3 



Z \ fcn 



7, 



O \ Km 



A- 



-2n B -3 / 3 / k x 2N 

r h + 2-; 2 U7 



Finally, the autocorrelation function for the Lorentz term is determined by, 

v L (k,k m ) = I ^E p B(g) Pi ;y:! 1) [fc 2 -2^+ g 2 +(fcx- g ) 2 

47r XoY l^-gl 2 

+4g(fex - g)(l - x 2 ) + 4<7 2 (1 - x 2 ) 2 ] 

which for the case at hand in the continuum limit leads to, 



V L (k,k n 



9 



[rm] 2 



PBO 

P70 



dxe v fe m/ 
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"in 



2(n B +3) 



dzz nB+2 e 



-1 



(l - 2zx + z 2 ) 2 [l + 2z 2 + (1 - 4z 2 )x 2 - 4zx 3 + 4z 



VI 



In the limit — <C 1 this can be approximated by, 



r L (k,k n 



Fmr 



PBO 

P70 
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A; 

.-(riB+|) 



-(n B +3) 



2(n B +3) 



/ ^B + 3 / A: 

\ Z \ /Cm 



+ — 2~^ B ' 2> 

15 Vfcn 



-2n B -3 / 3 , k x :^ 

r - + 2 ;2 U; 



(3.25) 



(3.26) 



(3.27) 



(3.28) 



It is interesting to note that the spectral functions only depend on the ratio k/k m . In figure [Q 
the correlation functions are shown. In figure [2] the spectrum determining the autocorrelation 
function of the anisotropic stress is shown for various values of the magnetic field strength and 
its spectral index. The approximate analytical solutions capture the general shape of the different 
power spectra, however, the amplitude is in all cases far below the one of the numerical solution. 
Only for wave numbers k much smaller than the one corresponding to the damping scale is the 
difference in the two curves not significant. 

4 Results 

Using the modified version of CMBEASY including the contributions due to a stochastic magnetic 
field present before decoupling whose spectrum is effectively cut-off at the diffusion scale using a 
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Figure 1: The spectra determining the autocorrelation functions of the magnetic anisotropic stress, 
energy density contrast and the Lorentz term are shown. For B = 10 nG and spectral index 
n-Q = —2.9 in each graph the numerical solution of the double integral, together with the numer- 
ical interpolation (spline) used in the code to calculate the CMB anisotropies and the analytical 
approximation in terms of incomplete Gamma functions is shown. The lower horizontal axis shows 
k, the upper one shows the ratio k/k m . In this example, the magnetic damping wave number is 
given by k m = 20 Mpc -1 (cf. equation (|3.15p ). 
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Figure 2: The spectra determining the autocorrelation function of the magnetic anisotropic stress 
V- KB are shown for different values of the magnetic field strength B and magnetic spectral index 
n-Q. In this case the numerical interpolation of the double integral appearing in equation (|3.24p was 
used. The magnetic damping wave number k m is given by equation (I3.15I) . 

gaussian window function (cf. equation (|3.12p ) the angular power spectra of the CMB temperature 
anisotropies and the polarization are calculated. In particular, for Pab > ^V B an d Vl the numerical 
interpolation as shown in figures [Hand [2] is used. This is much more accurate than the approximate 
solutions given in equations (|3.22[) . (|3.25p and (|3.28|) . but at the same time it is much faster than 
a full numerical integration for each step in A;-space. The CMB anisotropies have been calculated 
by using the initial conditions (|2.2ip as a whole including the curvature perturbation and the part 
involving the magnetic field. This corresponds to a complete correlation between the primordial 
curvature and the magnetic field [12j . Moreover, the primordial curvature perturbation is assumed 
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to be determined by the bestfit parameters of the ACDM model of WMAP7 [25]. In the numerical 
solution the values for Ab , 7Tb , L and £ were determined by the values of the square root of the 
corresponding power spectrum. In comparison to previous work [11] also the power spectrum for 
the Lorentz term is employed here which is more accurate since Ab and tt-q are correlated. This 
modification manifests itself in the interesting peak structure in the temperature anisotropies when 
comparing the cases with and without a magnetic field as discussed below. In figure [3] the angular 
power spectra of the temperature autocorrelation function, of the polarization and the temperature 
polarization cross correlation are shown for different values of the magnetic field parameters. These 
are compared with the angular power spectra obtained for the bestfit parameters of the ACDM 
model of WMAP7 [28], without a magnetic field. The magnetic field spectral index is chosen to 
be negative which is the case for magnetic fields generated during inflation (e.g. [23]) but not for 
those generated by a causal process, e.g. during the electroweak phase transition (e.g. [SI]), which 
requires tib to be an even integer and n-Q > 2 [35] . 
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Figure 3: The angular power spectra determining the autocorrelation functions of the temperature 
anisotropies of the CMB (Cj T ), the polarization of the CMB (Cf E ) and the temperature polar- 
ization cross correlation function (Cj ) for different values of the magnetic field strength while 
keeping fixed the magnetic field spectral index tt-b = —2.9. These have been calculated using the 
bestfit values of the six-parameter ACDM fit of WMAP7. For comparison, the pure ACDM model 
corresponding to B = has also been shown. 



Figures [3] and S] have been calculated using the bestfit values of the six parameter ACDM model 
of WMAP7 [28], that is in particular, Q b = 0.0445, fl A = 0.738, A^ = 2.38 x 10~ 9 , n s = 0.969 and 
the reionization optical depth r = 0.086. In figure [5] the contributions from the Sachs- Wolfe effect, 
the integrated Sachs- Wolfe effect and the Doppler term are shown for a model with a magnetic 
field and the bestfit ACDM model of WMAP7. 



4.1 The Sachs- Wolfe effect 

On large scales the temperature perturbation, written in the gauge invariant variables used here, 
is determined by \36\ [26] 
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Figure 4: The angular power spectra determining the autocorrelation functions of the temperature 
anisotropies of the CMB (Cj T ) for the magnetic field strength B = 5nG while varying the spectral 
index. These have been calculated using the bestfit values of the six-parameter ACDM fit of 
WMAP7. For comparison, the ACDM model with B = has also been shown. Left: The TT 
angular power spectrum for the whole range of multipoles. Middle: The first acoustic peak. Right: 
The second acoustic peak. 
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Figure 5: The different contributions to the total temperature angular power spectrum due to the 
Sachs- Wolfe term (SW), the integrated Sachs- Wolfe effect (ISW) and the Doppler term. The light 
solid line shows always the case in the presence of a magnetic field with B = 10 nG and n-Q = —2.9. 

where decoupling of the photons from the baryon fluid takes place instantaneously at the time of 
last scattering tls and tq corresponds to the present time. Moreover the Doppler term describing 
the relative motion between emitter and observer, namely Vbjn- 1 has been neglected since on large 
scales scales at last scattering it is negligible. This can also be appreciated in figure[5l The curvature 



14 



perturbation on slices of uniform total energy density in terms of the gauge-invariant variables used 
here is given by [16] 



c 



£ q a q o c 



I] Q 3(1 +W a )tta 

which on large scales is equivalent to equation (|2.22[) and also 



(4.30) 



c = $ 



3 1 + w 



(4.31) 



which during matter domination on large scales implies the standard result ( = |<3? (cf, e.g., [36J). 
The total curvature perturbation for the initial conditions (|2.2ip is given by 

D 7 A B 



4 Za 3(1 + w a )Q c 



(4.32) 



The photon density contrast A 7 is approximately constant on large scales (cf. equation (|2.1ip ) 
so that for the case of no magnetic field the curvature perturbation is conserved on large scales. 
However, in the presence of the magnetic it can be seen from equation ()4.32p that £ is approximately 
conserved on large scales during the radiation dominated epoch but not during matter domination 
where it behaves as £ = c± + C2T~ 2 , where c\ and C2 are constants. Assuming that last scattering 
takes place in the matter dominated era the curvature perturbation at tls is given by, 
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Therefore the gauge potential $ on large scales is given at the time of last scattering by 



$(tls) = jrCfa 



Q 7 (tj) ^ 7 (t ls ) 



A, 



(4.33) 



(4.34) 



where Tj refers to the time when the initial conditions are set deep inside the radiation dominated 
era. Clearly, in the case of no magnetic field C( r «) = C( t ls) and the standard result is recovered. Fi- 
nally, neglecting the integrated Sachs Wolfe effect (the integral in equation (J4.29J) ) the temperature 
perturbation is found to be 
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(4.35) 



Since f} 7 (Tj) is much larger than Q 7 (tls) it can be seen that the effect of the magnetic field 
is to reduce the temperature fluctuation on large scales. Thus the result of [TT] derived in the 
synchronous gauge (see also [37]) is recovered. This indicates that effectively the Sachs- Wolfe 
plateau is lower in the presence of the magnetic field than in the case without a magnetic field. 
This describes qualitatively the behaviour of the temperature fluctuation on large scales when 
initial conditions are set after neutrino decoupling. The case discussed here is different from the 
one discussed in [26J where the evolution since before the magnetic field generation which takes 
place before neutrino decoupling is taken into account. 
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4.2 Acoustic oscillations 



The effect of the magnetic field is strongly noticed in the acoustic oscillations of the baryon-photon 
fluid, as can be appreciated from the Sachs- Wolfe contribution in figure [5] (see also figures [3] and 
H]). These are determined by the photon energy density contrast A 7 . Using equations (I2.11h and 
(|2.15p its evolution is determined in the tight coupling limit by, 



A^,+ 
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where R h = -^ and c 2 b7 
solved by, [31] 
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tt j., : _ 1 is the sound speed of the baryon-photon fluid. In general this is 
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where in the case at hand F{r) = — ^ R , L+— A^ &. Furthermore, the soundhorizon is defined 
by t s (t) = L c s b 7 dr'. Neglecting the changes of i? b and $ due to the expansion of the universe 
|38j one arrives at a simple toy model to study the basic effects of the magnetic field. In this case 
equation (J4.36P reads, 



k 2 
[1 + i? b )A 7 + — A 7 
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which is solved by 



A 7 (r) = (A 7 (0) + L - 4(2 + i? b )$) cos(wt) + 
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to 



sinwr - L + 4(2 + R b )$, 



(4.38) 



(4.39) 



where oj = fc/y3(l + R\,). Here it can be seen that the contribution due to the magnetic field 
changes as well the amplitude of the oscillations in the photon energy density contrast as well 
as the zero point of the oscillations. Thus the effect of the magnetic field is similar to that of a 
change in the baryon energy density (cf., e.g., |31|). Similar to the case of the baryons here as well 
an alternating peak structure is found, that is whereas the odd-numbered peaks are lower than 
in the case without a magnetic field, the even-numbered peaks are larger. In figure [6] numerical 
solutions for —?- + \P — $ for several wave numbers are shown. The general behaviour can be 
traced back to the effect the magnetic field contribution has on the oscillations in the photon 
energy density contrast. Going back to equation (J4.36P and bearing in mind that L is positive 
in the model at hand it can be seen that the magnetic field effectively augments the pressure 
perturbation of the baryon-photon fluid. At last scattering there are mainly two competing effects 
coming from the gravitational infall determined by the gravitational potential ^ and the restoring 
force due to pressure determined in the case at hand by c 2 b7 A 7 . The odd-numbered acoustic peaks 



correspond to the compressional phase of the oscillations [31J. Thus increasing the pressure while 
not changing the gravitational potential significantly results in a larger restoring force and so less 
effect of the gravitational infall. This reduces the amplitude of odd-numbered peaks in the models 
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Figure 6: The numerical solution of the evolution of the effective temperature perturbation is shown 
for different wave numbers. The WMAP7 bestfit solution is shown in comparison with the solution 
in presence of a stochastic magnetic field with B = 20 nG and magnetic spectral index n-Q = —2.9. 
Indicated is the time of last scattering, which in this case is tls = 285 Mpc. Left: The solutions 
for the wave numbers k = 0.0140 Mpc -1 and k = 0.0172 Mpc -1 have a maximum at the time of 
last scattering. These wave numbers correspond to the region of the first acoustic peak, whereas 
the former corresponds to I = 200, the latter corresponds to £ = 244. Right: The solutions for 
the wave numbers k = 0.0364 Mpc -1 and k = 0.0396 Mpc -1 have a minimum at the time of last 
scattering. These wave numbers correspond to the region of the second acoustic peak, whereas the 
former corresponds to £ = 523, the latter corresponds to £ = 569. 

with a magnetic field (cf. figured]). On the contrary, the even numbered peaks correspond to the 
expansion phase inside the gravitational well [31]. Since the pressure is larger in models with a 
magnetic field the resulting even-numbered peaks are larger than in the model without a magnetic 
field. This situation is reflected in figure [6] for the second peak which corresponds to a minimum 
in A 7 at the time of last scattering. As can be seen in figure [6] the minimum in the presence of a 
magnetic field is deeper than in the case without. Thus leading to a larger amplitude of the second 
acoustic peak in the presence of a magnetic field. This alternating peak structure is clearly seen 
in the numerical solutions, cf. figure [3]l5j Moreover, these characteristics are present as well in 
the extrema of the toy model (|4.39[) . Assuming that the initial value of the effective temperature 
perturbation = -j + ^ — <J>, once it enters the horizon, is of the order of its Sachs- Wolfe amplitude 
the extrema of O can be estimated. At wtls = nir, corresponding to a maximum, it is found that 
@max = 3(1 + 6.Rb)<J?(TLs) ~~ 2 an< ^ a ^ ^^s = 2n7r, corresponding to a minimum Q m in = — ^(tls)- 
So that in general the maxima are expected to be smaller in the presence of a magnetic field which 
results in comparatively smaller values of the odd-numbered acoustic peaks. 
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4.3 The linear matter power spectrum 



During matter domination the total matter perturbation A m = i? c A c + i?b^b) where Ri determines 
the fractional energy density with respect to the total matter energy density, is determined by, 



A m + HA r 
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In deriving equation (J4.4Q[) it has been used that \& ~ — <J> and that during matter domination on 



small scales, <& 



2M£k 2 



which yields to <£ 



3H 2 

Ik? 



(A m + S7 7 Ab). Thus on very small scales, the 



contribution from the magnetic field on the righthandside of equation (|4.40p dominates. In this 
limit there is a constant solution of equation (|4.40p . namely, A m oc k 2 L. This implies that for large 
values of k the matter power spectrum scales as the corresponding power spectrum of the Lorentz 
term, Va^ oc k 4 Vi, [33]. The total matter power spectrum VA m is shown in figure [7] which clearly 
shows the influence of the magnetic field on small scales. The resulting linear matter power spectra 
are similar to what is presented in |12| . However, the general shape differs from that reported 
in |33| which only considers the compensated magnetic mode without the contribution from the 
primordial curvature mode and uses a different form for the power spectra related to the magnetic 
contribution. The normalization of the matter power spectrum as is defined to be the variance of 
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Figure 7: The total matter power spectrum for different spectral index and amplitude of the 
magnetic field. The data points included show the estimated power spectrum from the APM 
angular galaxy catalogue [40] (see also [IT]). 



a mass fluctuation inside a sphere of radius R = 8h 1 Mpc, that is (see for example, |39|) 
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(4.41) 



where ji(x) = (sinx — xcosx)/x 2 . Calculating as for the different values of the spectral index 
and strength of the magnetic field used in the examples above it is found that for n B = —2.9 ag 



18 



decreases with the strength of the magnetic field, so that for B = 10 nG it is found that as = 0.704 
whereas for the ACDM model without magnetic field as = 0.802. On the other hand varying the 
spectral index and keeping the magnetic field strength fixed at B = 5 nG the largest value is found 
for n B = —1.5 that is eg = 0.802 and the smallest for n B = —2.9 with as = 0.77. 



4.4 Uncorrelated curvature and compensated magnetic modes 

Before closing this section we will comment on the case in which the primordial curvature mode, 
which is assumed to not include any magnetic contributions, and the compensated magnetic mode 
are uncorrelated. The passive mode is not considered here, since it was also not considered explicitly 
in the previous subsections where the curvature perturbation was assumed to be determined by the 
bestfit values of the six parameter ACDM model of WMAP7 [28] . 

The compensated magnetic mode corresponds to zero inital curvature [TJj and we follow the 
treatment of correlated isocurvature initial conditions (see for example, |42|). In this case the initial 
conditions (cf. equations (|2.2ip ) with vanishing total curvature perturbation (£ = 0) have two 
different contributions, one proportional to the magnetic energy density Ab and one proportional 
to the magnetic anisotropic stress tt b . Moreover, the Lorentz term can also be written in terms 
of Ab and 7Tb, (cf. equation (|3.8p ). Thus the initial conditions as well as the evolution equations 
only contain terms proportional to the magnetic energy density or the magnetic anisotropic stress, 
respectively. Therefore we define the total brightness function 0^ due to the compensated magnetic 
mode by 



Qt(k) = G* B (k)A B (k) + G^{k)^ B {k), 



where G^ B (k) is the transfer function calculated for Ab 



(4.42) 



1 and 7Tb = and G^ B (k) is the transfer 
function obtained for Ab = and 7Tb = 1. Moreover, Ab denotes the Gaussian random variable 
corresponding to the magnetic energy density and n B the Gaussian random variable corresponding 
to the magnetic anisotropic stress. The autocorrelation functions of these variables are determined 
by equations (|3.19|) . (|3.2ip and (|3.24|> . In addition, the cross correlation function of Ab and % B is 
required 



where 



^Ab?tb (fc) 



dxe 



. - ->, 2vr 2 

(A B (k)fr B (k')) = -p-^A BWB (fc)%, 



[rm] 



+3M 2 



PESO 

P70 



±Y CnB+ Vte) s 

A',, 



n-a—1 



'l-2zx + z 2 ) 2 (-l + z 2 + 



zx 



dzz nB+2 e 2 



1 + 3z 2 )x 2 + 3zx 3 



(4.43) 



(4.44) 



which is shown in figure [8] for particular values of the magnetic field parameters. The total angular 
power spectrum determining the temperature autocorrelation function, Cf is given by, 
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Figure 8: The spectrum determining the crosscorrelation function of the magnetic energy density 
and anisotropic stress is shown for B = 10 nG and spectral index n-Q = —2.9 . The numerical 
solution of the double integral, together with the numerical interpolation (spline) used in the code 
to calculate the CMB anisotropies is plotted. The lower horizontal axis shows k, the upper one 
shows the ratio k/k m . In this example, the magnetic damping wave number is given by k m = 20 
Mpc" 1 (cf. equation ([37151) ). 



and a similar expression for the angular power spectrum determining the autocorrelation function 
of the E-mode, Cf . The temperature-polarization cross correlation angular power spectrum is 
given by, 
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(4.46) 



where H e B and H^ B denote the polarization transfer functions for Ab = 1, 7Tb = and Ab = 0, 
7Tb = 1, respectively. The angular power spectra due to the compensated magnetic mode are shown 
in figure [9j Assuming that there is no correlation between the magnetic mode and the primordial 
curvature fluctuation the total angular power spectrum is given by the sum of the angular power 
spectra. Taking into account the order of magnitude of the bestfit WMAP7 ACDM model (cf. 
figure [3]) and the compensated magnetic mode, as shown in figure [9] for B = 10 nG and n-g = —2.9, 
it is seen that the contribution of the magnetic field to the final angular power spectrum is much 
smaller than when treating the curvature perturbation and the magnetic mode as fully correlated 
as in the previous subsections (cf . figure [3]) . 

Finally, the total linear matter perturbation due to the compensated magnetic mode is given 
by 
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Figure 9: The angular power spectra determining the autocorrelation functions of the temperature 
anisotropies of the CMB (Cj T ) (cf. equation (|4.45[> ). the polarization of the CMB (Cf E ) and the 
temperature polarization cross correlation function (Cj ) (cf. equation (|4.46p ) for the compensated 
magnetic mode for B = 10 nG and n-Q = —2.9. These have been calculated using the bestfit values 
of the six parameter ACDM fit of WMAP7. 

where A^Ab an d Aj i7rB denote the final density perturbations proportional to Ab and 7Tb, respec- 
tively, and i denotes b or c. In the former case only the contribution due to the magnetic energy 
density is taken into account in the initial conditions and the evolution equations and in the lat- 
ter only the one due to the magnetic anisotropic stress. Therefore, the total linear matter power 
spectrum due to the compensated magnetic mode is found to be 

"" \RlU h + 2R b R c V + R 2 C U C 



PAjk) 



A- 3 



where 



Ui 
V 



Pa b A? >Ab + 2V Ab7TB Ai,A B A^ B + 7\ B A^ i = b, c 

"PAbA^AbAcAb + - PAb7tb (Ab,A B Ac,7r B + Ab l7rB A Ci A B ) + "^ttb Ab,7r B A 



C,7TB! 



(4.48) 



(4.49) 



where the definitions A^Ab = AbAj ) Ab=i,7t b =o and Aj j7I - B = vrBAi j7rB= i i A B =o were used. The total 
matter power spectrum due to the compensated magnetic mode is shown in figure [10] (left) for 
different values of the magnetic field parameters which is similar to the result of |33j, but not 
identical due to the gaussian window function used here in the calculation of the different magnetic 
spectra. Assuming that the standard adiabatic mode, due to a primordial curvature perturbation, 
and the magnetic mode are uncorrelated the total matter power spectrum for the WMAP7 ACDM 
model is given by the sum of the matter power spectra which is shown in figure [10] (right). In 
this case the resulting value of a$ including the contributions from the pure adiabatic mode of the 
bestfit WMAP7 ACDM model and the compensated magnetic mode is always larger than in the 
case with no magnetic field. 

5 Conclusions 

The CMB anisotropies in the presence of a stochastic magnetic field have been calculated using the 
gauge invariant formalism. Initial conditions have been found by imposing adiabaticity between 
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Figure 10: Left: The total linear matter power spectrum due to the compensated magnetic mode. 
For comparison the total linear matter power spectrum for the WMAP7 bestfit model is included 
which corresponds to a pure adiabatic mode. Right: The total linear matter power spectrum 
due to the adiabatic mode and the compensated magnetic mode assuming that these are uncorre- 
cted. Both graphs have been calculated using the bestfit values of the six parameter ACDM fit of 
WMAP7. 

all matter /radiation components except the magnetic contribution. This ensures that in the case 
without a magnetic field the standard adiabatic mode is recovered. The numerical solutions were 
calculated with a modified version of the CMBEASY programme [9]. The magnetic field is assumed 
to be a gaussian field whose spectrum contains a gaussian window function effectively cutting off 
the spectrum at a wave number corresponding to the magnetic damping scale. The magnetic field 
contributions in the evolution equations require the corresponding power spectra of the magnetic 
energy density, the anisotropic stress, the Lorentz term and the cross correlation between the 
magnetic energy density and magnetic stress. These involve convolution integrals which only in the 
limit of scales much larger than the damping scale can be well approximated in terms of special 
functions. Therefore the numerical solutions have been obtained using a numerical interpolation 
for these integrals. A strategy similar to the one used in the different Boltzmann solver codes and 
in particular CMBEASY to calculate the angular power spectrum of the CMB anisotropics. 

Firstly, based on previous work |llj . the initial conditions are used as a whole, including the 
curvature mode and the part of the magnetic field. Effectively, this corresponds to a complete 
correlation between the primordial curvature perturbation and the magnetic field. The curvature 
mode is assumed to be determined by the bestfit values of the six parameter ACDM model of 
WMAP7 [28j. The magnetic field variables are determined by the square root of the corresponding 
power spectra. The difference to the numerical solutions presented in |llj is that here the power 
spectrum for the Lorentz term was used which is more accurate since the magnetic energy density 
and anisotropic stress are correlated. Moreover, no approximations have been used for the power 
spectra related to the magnetic field contributions. Having obtained the numerical solutions for 
different values of the magnetic field amplitude and spectral index the results were compared with 
qualitative estimates on very large scales and on small scales. In particular the Sachs- Wolfe effect 
was calculated which was found to be smaller in the presence of a magnetic field which is also what 
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is observed in the numerical solutions. Moreover the acoustic oscillations in the photon-baryon 
fluid were described in a simple model and the alternating peak structure recovered which is seen 
in the numerical solution for the type of magnetic field spectra used here. 

Secondly, the case of uncorrelated curvature and compensated magnetic modes is considered. 
The compensated magnetic mode is considered to be a correlated isocurvature mode consisting of 
two modes, one proportional to the magnetic energy density and one proportional to the magnetic 
anisotropic stress. The final angular power spectra determining the temperature and polarization 
auto- and cross-correlations are calculated taking into account the cross correlation between the 
magnetic energy density and anisotropic stress. Assuming no correlation of the compensated mag- 
netic mode with the standard adiabatic mode due to a primordial curvature fluctuation the overall 
effect of the magnetic field contribution is much smaller than for the same magnetic field values in 
the fully correlated case. Moreover, the total angular power spectra resulting from the standard 
adiabatic mode and the compensated magnetic mode are always larger in the case including the 
magnetic field. 

Finally, the linear matter power spectrum was calculated in both cases. At small scales it is 
observed that it is dominated by the magnetic field contribution which is also seen in the second 
order differential equation determining the total matter density contrast derived in this limit. The 
total linear matter power spectrum is found to have distinct shapes in the two different cases 
assuming correlated or uncorrelated curvature and magnetic modes, respectively. Moreover, it is 
also noted that the normalization of the matter density fluctuation og depends on the magnetic field 
parameters. This provides yet another possibility of putting limits on the magnetic field parameters 

m. 
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7 Appendix: Transformation of initial conditions to the synchronous 
gauge 

The initial conditions (cf. equations (|2.2ip ) are given in terms of the gauge- invariant variables. 
The corresponding initial conditions in the synchronous gauge are found by using the following. In 
general the perturbed metric is given by |16t [i~5] , 

ds 2 = a 2 {r) [-(1 + 2A)dr 2 - 2B l drdx i + {5 lj + 2H ij )dx l dx j ] (7.50) 

where the expansion in harmonic functions Y(k, x), as defined in section 3, for scalar perturbations 
leads to 

Bi = BYi H i3 = H L Y6 l:J + tf T l%. (7.51) 

In the synchronous gauge A = 0, B = 0. In this case, it is common to write the metric perturbation 
as EH, 



hij(x,T) = I d 3 ke ik -* 



kikjh(k,T) + (kikj - -%)6r/(A;,r) 
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k = j, (7-52) 



which implies rj = — (H-^ + \H<?) and h = 6#l- The gauge invariant Bardeen potentials ^ and <3? 
are given, respectively, by |16l I15j . Vf = A — 1-Lk~ l a — k~ x & and <!> = Hi, + |#t — Tik~ 1 a, where 
a = A; _1 i7x — B. Moreover, the gauge-invariant velocity is defined as V = v — /c -1 £Tt- Thus the 
velocity in the synchronous gauge is found to be, 

v s = V - -A- I dT'a(T')^>(T r ). (7.53) 

a(r) J 

Therefore together with equation (12. 4D the initial conditions (I2.21D in the synchronous gauge are 
given by, 

«? = ^ = ^ = ^ = -^(ab + L + ^vtb)=-^ 7 Ab, 

7 b 8 V 3 / 8 46 

o f/L / . ^2 $Xv + 1 \ O-v . H-v ^b 

^ = -f ( AB + L -3^ B J=-f AB + rflP 

S c 5 = 

fL _ 4C ^ 4vr B - 60„A B 



n v 15 + 4n„ 2Q„ 15 + 4ft 



z/ 



V = -C + ^(-Ab + L + ^7Tb) =-C, (7.54) 

where v s = v s /x. There are two integration constants one in the calculation of v s and one in the 
calculation of h. These have been chosen such that v^ = and there is no constant mode in h. 
The final expression on the far right of each equation is obtained using L = Ab — q^b- For ttq only 
the expression using this relation for L is given. At this order of x, h = 0. These are the initial 
conditions to lowest order in x in the synchronous gauge as given in |18^ \TT\ [T2] . 
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